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Abstract 



Quantum corrections to the properties of a homogeneous interacting Bose gas at 
zero temperature can be calculated as a low-density expansion in powers of \/pa^, 
^ . where p is the number density and a is the S-wave scattering length. We calculate 

the ground state energy density to second order in \/pa^. The coefficient of the pa^ 
correction has a logarithmic term that was calculated in 1959. We present the first 
^SJ . calculation of the constant under the logarithm. The constant depends not only on a, 

but also on an extra parameter that describes the low energy 3 — > 3 scattering of the 
bosons. In the case of alkali atoms, we argue that the second order quantum correction 
is dominated by the logarithmic term, where the argument of the logarithm is paiy, 
and iy is the length scale set by the van der Waals potential. 



1 Introduction 



The successful achievement of Bose-Einstein condensation of atomic gases in magnetic traps 
|1|, ^ 1^ has created a revival of interest in Bose gases. While a qualitative description of the 
condensation can be obtained using mean field methods [^], a more quantitative treatment 
requires including corrections from quantum fluctuations around the mean-field. The relative 
magnitude of these corrections grows as the square root of the number density of the atoms. 
They will therefore become more important as higher condensate densities are achieved. 

In order to develop a deeper understanding of the quantum fluctuations, it is worthwhile 
to go back to the simpler problem of a homogeneous gas of interacting bosons at zero 
temperature. This problem was studied intensively in the 1950's p, 0. The properties 
of the system can be calculated as an expansion in powers of a/po^, where p is the number 
density of atoms and a is their S-wave scattering length. For example, the expansion for the 
energy density has the form 
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The coefficient of y/pa-^ was first obtained by Lee and Yang for a hard sphere gas |^. The 
pa^ correction is the first term in the expansion that is sensitive to atomic parameters other 
than the scattering length. The coefficient of ln(p) in the pa^ correction was calculated by 
Wu, by Hugenholtz and Pines, and by Sawada in 1959 0. 

In this paper, we present the first calculation of the constant C under the logarithm in 
(|I|). We express the constant in terms of a coupling constant g^ that is defined in terms of 
the low-energy behavior of the amplitude for the 3^3 scattering of atoms in the vacuum. 
The scattering length a and the coupling constant g^ are the only atomic physics parameters 
that are needed to calculate the energy density to second order in the quantum corrections. 
Since g^ is difficult to measure experimentally or calculate theoretically, it must be treated 
as a phenomenological parameter. The dependence on this undetermined parameter creates 
a large uncertainty in the second order quantum correction, except in cases where the cor- 
rection is dominated by the logarithmic term. We argue that this will typically be the case 
for alkali atoms, provided we take the argument of the logarithm to be paiy, where iy is 
the length scale set by the van der Waals potential. In this case our result reduces to 
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If the logarithm is large compared to 1, this should give an accurate estimate of the second 
order quantum corrections to the energy density. 

We begin in Section |^ by formulating the problem of calculating the energy density as 
a quantum field theory problem. In Section ^, we set up a perturbative framework and use 
it to calculate the energy density to second order in the quantum corrections. In Section ^ 
we calculate the second order quantum correction to the T-matrix element for the 3-^3 
scattering of atoms in the vacuum. We define the coupling constant g^ and show that the 
renormalization of g^ is necessary to remove a logarithmic ultraviolet divergence from the T- 
matrix element. In Section ^ we complete the calculation of the energy density by including 



the contribution from g^ and its renormalization. We then discuss the case of alkah atoms, 
and argue that the dependence on g^ can be ehminated in favor of a logarithmic dependence 
on the length scale set by the van der Waals interaction. The two-loop Feynman diagrams 
that contribute to the energy density are calculated in the Appendix. 

2 Field Theory Formulation 

We are interested in computing the ground state energy density £^ of a homogeneous Bose 
gas as a function of its density p. Our starting point is a local quantum field theory that 
describes atoms with momenta much lower than the inverse of the range of the interatomic 
potential, which is several A in the case of alkali atoms. At such low energies, the interactions 
appear pointlike on the scale of the de Broglie wavelengths of the atoms. The many-body 
quantum mechanics of the atoms can therefore be formulated in terms of a quantum field 
theory whose hamiltonian density is a local function of the field: 

n = |;^W^-W + ^^(^V)' + ^^3(^V)' + ■••• (3) 

For simplicity, we have assumed that the atoms have only one spin state so that they can be 
represented by a single complex field ip{'x,t). The {ip'^ip)'^ term represents 2 — i> 2 scattering 
through an S-wave interaction with scattering length a given by 

9 = • (4) 

m 

This coupling constant contains all the information about atomic interactions that is re- 
quired to calculate the first order quantum corrections to the properties of a sufficiently cold 
and dilute Bose gas. We follow the effective-field-theory philosophy [^ of including in the 
hamiltonian all possible local terms that are consistent with the symmetries, which include 
Galilean invariance and the phase symmetry ip — >■ e*"?/'. The term (t/'^"?/^)^ in (^ allows 3 — > 3 
scattering through a pointlike interaction. The . . .'s in (^) include all possible terms that are 
higher order in the derivatives or in the number of fields. In principle, the coefficients of these 
terms can be calculated from the n-body potentials that describe interatomic interactions. 
In the absence of such calculations, they can be taken as phenomenological parameters. 
The effective-field-theory philosophy is based on the assumption that there is a systematic 
expansion in powers of the momentum. The relative importance of the interactions terms 
in d^) at a low momentum scale p is then given by the dimensionless combination of the 
coupling constant and p, which is mgp and rngsp"^ for the {tp'^ipy term and the {ip^ip)^ term, 
respectively. If p is sufficiently small, the effects of the g^ term will be much smaller than 
those of the g term. Terms with more derivatives or with higher powers of ip give effects 
that are suppressed by even more powers of p. Of the infinitely many terms in (|^), there 
are only a finite number that contribute at any given order in p. In the case of the energy 
density, the appropriate momentum scales are p^^^ and (paY^'^, so the momentum expansion 
becomes an expansion in powers of the density. We will find that the scattering length g 
and the coupling constant g^ are the only atomic physics parameters that contribute to the 
energy density through third order in p [R|. 



At a given order in p, only a finite number of terms contribute. By tuning tlie coefficients 
of tliese terms, one can describe n ^ n scattering of atoms in tlie vacuum witli wliatever 
accuracy is desired. 

Tlie phase symmetry tp — > e*°^ of the hamiltonian implies the conservation of the number 
of atoms. The number density operator is A/" = ip^ip. A homogeneous Bose gas can be 
described by a field theory with hamiltonian density 7i — ^M, where /i is the chemical 
potential. The energy density £{^) and the number density p(yu) in the ground state of this 
field theory are 

£{^A = {n),, (5) 

p(/^) = (A/")., (6) 

where (■ ■ ■)^ denotes the expectation value in the ground state. By eliminating /x from the 
two equations (^) and (^, we obtain £^ as a function of p. It is simpler in perturbative 
calculations to first calculate the free energy density J^(/i) of the ground state: 

T{p)^{H-^lM),. (7) 

After inverting (|^) to obtain yU as a function of p, one can obtain the energy density from 

S = r + pp. (8) 

The partition function for the field theory with hamiltonian density Ti — pM can be 
expressed as a functional integral: 

Z = f ViP^ Vtfj exp {z S[^]} , (9) 

where the action S[ip] is given by 

S[^] = Jdt J d\ Lt U +Y_+f\^_].^g (^t^)2 _ ^^3 (^t^)3 _ \ (10) 

We have set /i = 1 in the action. Dimensional analysis can be used to reinsert the factors of 
h at the end of the calculation. The . . .'s in ([l0|) represents all possible terms with higher 
powers of ip or more factors of V. They also include counterterms that are needed to cancel 
ultraviolet divergences associated with the parameters /i, g, and g^. The free energy density 
jF(/i) is related to the partition function by 

Z = exp {-iVT T{p)} , (11) 

where VT is the spacetime volume. The ground state expectation value of an operator (O) 
can be expressed as a functional integral: 

(O), = ^Jv^^V^OexpitSm. (12) 

By differentiating the logarithm of both (^ and (|TT]) with respect to p and using (|I2]), we 
obtain the relation 

dp 



Differentiating (j^) with respect to p and using (|T^), we obtain 

^(p)=Mp)- (14) 

The simplest way to calculate the energy density is to first calculate J-'{fi), use (|T3D to get 



p(/x), invert to get /i(p), and then integrate (|Tj) to get £{p) 



It is convenient to parameterize the quantum field ^(r, t) in terms of two real-valued 

quantum fields ^ and r] that describe quantum fluctuations around an arbitrary constant 

background v. 

i( +\ , ^(r,^)+^^(r,t) 
^(r,t) = V + -= . (15) 

After inserting the field parameterization ( pTsD into the action (|iy), it can be expanded in 
powers of the quantum fields ^ and rj. By separating the action into a free part and an 
interaction part, we can express the thermodynamic functions as diagrammatic expansions. 
The free energy density T is the sum of all connected vacuum diagrams, which are Feynman 
diagrams with no external legs. This sum is independent of the arbitrary background v. It 
is convenient to define the thermodynamic potential r2(/i, f), which is the sum of all one- 
particle-irreducible vacuum diagrams. The thermodynamic potential, which depends on v, 
contains the information required to determine all of the thermodynamic functions. The free 
energy J-'{p) can be obtained by evaluating Q{fi, v) at a particular value of v given by the 
tadpole condition 

v{p) = {^),. (16) 

For this value of v, those diagrams that can be disconnected by cutting a single line vanish. 
Thus the sum of connected vacuum diagram reduces to the sum of one-particle-irreducible 
vacuum diagrams and we have 

^(^) = fi(/i,TJ(/i)). (17) 

Using (0), the tadpole condition ([T6|) reduces to (,^)^ = {r])^ = 0. The phase of the field ip 
can be chosen so that (r^)^ is automatically 0. The condition (^)^ = can be conveniently 
expressed in terms of the thermodynamic potential itself: 

^{pMf^)) = 0. (18) 

Differentiating both sides of (|l3) with respect to p and using (|^), we obtain 

Comparing with (|l^), we find that the number density can be expressed as 

P(P) = -^(/i,tJ(/i)). (20) 



3 Ground State Energy Density 

In this section, we calculate the ground state energy for a homogeneous Bose gas to second 
order in the quantum corrections. We first set up a perturbative framework for carrying 
out calculations in the presence of a nonzero chemical potential. We use the framework to 
calculate the energy density to second order in the quantum corrections. We then carry out 
the renormalizations of ^ and g that are necessary to remove power ultraviolet divergences 
from the energy density. 

3.1 Perturbative framework 



We can describe a Bose gas with nonzero density p by the action ([T0|) with an appropriately 
chosen value of the chemical potential. For simplicity, we set g^ = and omit all terms in (|TUp 
that are higher order in ip or V. We ignore for the moment the counterterms associated 
with renormalization, so the parameters /i and g should be regarded as bare parameters. 
Inserting the field parameterization ([T5|) into the action and expanding in powers of ^ and 
7], the action becomes 

S|^l = / * / d'xLv' - igi." + ^f + i (,,{ - {,)) + i-f (V^ - 2mgv' + x) ( 

where / = |^/ and 

X = 2m(^p- ^gv^^ . (22) 

To organize the quantum corrections into a loop expansion, we separate the terms in the 
action (^) that depend on ^ and t] into a free part and an interaction part: 

S[ij] = S[v] + 5frec[^,r/] + Si^,[v,^,r]]. (23) 

The free part of the action consists of the terms quadratic in C, and rj: 



^f,ee[e,^] = Jdtjd^xl^[r]^-^r] 



+ 4^e(V^ - 2mgv' + X)^ + ^r^iV' + X)r?| . (24) 

The Fourier transform of the propagator for the fields ^ and r/ is a 2 x 2 matrix: 

n/ i, ^_ i /(P-X)/2m -iu \ , . 

^^'^' "' ''^~ u^- e^{k, v)+ie[ 100 (P + 2mgv^ - X)/2m ) ' ^^^^ 

where k is the wavevector, uo is the frequency, and 

e^{k, v) = -i^(A;2 + 2mgv^ - X){k^ - X) . (26) 



The diagonal elements of the propagator matrix (pq) are represented by solid lines for ^ 
and dashed lines for rj, as illustrated in Figs. 1(a) and 1(b). The off-diagonal elements are 
represented by a line that is half solid and half dashed, as in Fig. 1(c). All the remaining 
terms in the action f^Tf) are treated as interactions: 



The term proportional to ^ is represented by a dot at which a solid line terminates as 
illustrated in Fig. 1(d). The 3 and 4-point couplings are represented by points that connect 
three and four lines, respectively, as in Figs. l(e)-l(i). 

It is possible to diagonalize the propagator matrix ( |25D by applying a Bogoliubov trans- 
formation to the fields C, and r]. However, such a transformation makes the interaction terms 
in the action significantly more complicated and increases the number of diagrams that 
contribute to most quantities. For explicit calculations, it is more economical to minimize 
the number of diagrams. We therefore prefer to use a propagator matrix with off-diagonal 
elements. This perturbative framework was recently used by Haugset, Haugerud, and Ravn- 
dal [§ to reproduce the -y/po^ correction in the expression (P for the energy density. 

3.2 Free energy density 

If the n-loop contribution to the thermodynamic potential Q is denoted by Qn{f^, v), the loop 
expansion for the free energy density ( piTf ) is 

^(/i) = (]o(/i,^) + fii(/i,^) + ^2{l^,v) + ■■■ , (28) 

where v is the condensate, which satisfies (ITSf): 



^(P,U) + ^(/^,^) + ■■■ = 0. (29) 

The loop expansion (^) does not coincide with the expansion in the order of quantum 
corrections because of its dependence on v. To obtain the quantum expansion, we must 
expand the condensate v around its classical value wq, which satisfies 

^(/U,tJo) = 0. (30) 

By expanding ( pUD in powers oiv — tJo, and solving for U, we obtain the quantum expansion 
for the condensate: 

v{li) = vo{ft) + vi{fi) + V2{fi) + ••• , (31) 

where U„ is the n-th order quantum correction. For example, the first-order quantum cor- 
rection is 



Inserting (|3T|) into ( pS]) and expanding in powers of fi, f2, . • ., we obtain the quantum 
expansion for the free energy. Keeping only terms through second order, we have 



The mean-field contribution to fi(/i, v) is given by the terms in (21) that are independent 
of C, and 7]-. 

noif^,v) = -i^v^ + ^gv\ (34) 

One solution to ( PDD is vq = 0, but it is a local maximum of Qq and therefore represents an 
unstable configuration. The stable solution is 

vl = 2fi/g. (35) 

The tree-level contribution to J^(/i) given in ( ^3]) is 



^o{lJ^,Vo) = . (36) 

9 

The dispersion relation (p6|) simplifies significantly at the point v = Vq, because X = at 
that point. It reduces to the Bogoliubov dispersion relation: 



e{k) = s{k,vo) = , (37) 

2m 

where A^ = 4m/i. The tadpole interaction in (pT]) also vanishes when v = Vq. 

Using the free part of the action (p^ , we can obtain the one-loop contribution to fi(/i, v): 

Q,{^,v)='-J-0^J^\ndetD{u;,k,v), (38) 



where D(uj, k,v) is given in (^Sf ). By integrating over u, we obtain 

1 r d^k 
^i(/^,^) = 2 7 7^^(^'^)' (39) 

where e{k,v) is given by (|26|). The one-loop contribution to the free energy is 

fii(/i,Uo) = -—Io,-ii^mfx) , (40) 

Am 

where /o,-i is a function of A^ = 4m,fi defined by the integral 

/^„(A2) = f ^\ (pT . (41) 

Differentiating the expression (p9D with respect to f and evaluating at f = IJo, we obtain 

^(^«„) = f(3/u + /-w). (42) 



The first order quantum correction to the condensate, which is given by (|32D is 

gvo 



Vl 



16/i 



(3/: 



1,1 



^-1,-1 



(43) 



The two-loop contribution to the thermodynamic potential is obtained from the vacuum 
diagrams shown in Fig. ^. The contributions from the individual diagrams are given in 
Appendix |A.5| . The sum of the diagrams gives 



^2{fJ',Vo] 



8 64 



3/^ _i + 2/_i,_i/i,i + 3/^ 



where 



J — 6 Jo,o,i ~ ^-1,-1,1 ~ 3 Ji^i^i — 2 J_i 0,0 • 
The integrals Ji^rn,n are functions of A^ = 4myU defined by 



■Jl,m,n[,^'- J 



d^p f d?q [p^ /2me{j))\ [g^/2me(g)]™' [r^/2m£(r)]' 



{2Tif J (27r)3 



2m[e{p) + e{q) + e{r)\ 



(44) 
(45) 

(46) 



where ^ = |p + q|- 

Inserting (|36D , (^0|), (l4^), (143|) , and (|^) into (|33|), we obtain the complete expression for 
the free energy density to second order in the quantum corrections: 



Tr( \ ^^ < ^ T (A ^ , ^3^^ T < 9 \j2 oT T 

^{N = + ^^0,-1 4m/i + ^— J + — /_! 1 - 2/_i,_i/i,i 

g Am 8 62 <- ' 



or2 

-J -'1.1 



(47) 



The free energy ( P7[ ) depends on fi through the explicit factors of fi and through the momen- 
tum scale of the integrals which is A^ = 4m/x. We have made the argument of the integral 
explicit for the /o,-i term in (W^). 



3.3 Energy density 

To calculate the energy density, we use (p!3| ) to get p{fi), invert that relation to get fi{p), and 



integrate using (]T^) to get £{p)- Differentiating (^) and using (|13]), the number density is 



p(/i) = ^ - i/i,i(4m/i) - ^ [J + 4m/iJ' + (/_i,_i - /i,i)/o,i + (/-i,-i + 3/i,i)/2,3] , (48) 
g z » 

where J' = dJ/dA"^ and all the integrals are functions of A^ = Amfi. We have made the 
argument of the integral explicit for the Ji,i term. We have used the identity ( |A^) to 
differentiate the integrals Im,n- Inverting (^Sf ) to get /x as a function of p, we obtain 

1 1 . ,_ ^ mg"^ 



M = 1^9 P + -^ghA^mgp) 



16 



[J + 2mgpJ' + (/_i,_i - /i,i)/o,i + (/-i,-i + /i,i)/2,3] , 

(49) 
we can write 



where the integrals are now functions of A^ = 2mgp. Using the identity 
the expression as a total derivative: 



Kp) 



d r 

d~p\A 



1 
4m 



9P^ + -;-Io,-ii'^mgp) + 



mg^P j,9_ 
16 +32 



-^-1,-1 ~ 2J_i,_iJi,i — III 



(50) 



We can now read off the energy density using ( n ) : 



8{p) = Sq + -gp^ + —Io,-i{2mgp) + 



16 32 L -I'-i ^' ^'^ I'l 



(51) 



where So is an integration constant and all the integrals are functions of A^ = 2mgp. It is 
convenient to choose the integration constant Sq so that the energy of the vacuum is zero: 
^(0) =0. 



3.4 Renormalization of /i and g 

Our result (^) for the energy density can be generalized to an arbitrary number of spatial 
dimensions D simply by replacing the integration measure (i^p/(27r)^ in (^Tj) and p6| ) by 



d^p/{27i)^. The integrals Im,n and Ji^rn,n in (|5lD are ultraviolet divergent for any positive 
number of dimensions D. If we impose a momentum cutoff Ajjv, then /o,-i diverges like 
A^y^, while /-i,-i and Iii diverge like A^y. The integrals Ji^rn,n contain subintegrals that 
diverge like A^y^, and, if D > 1, they also have an overall divergence that scales like A^^~^. 
There are cancellations among the Ji^rn,n integrals that reduce the overall divergence to A^^~^ 
for D > 3 and InAf/y for D = 3. 

The divergences can be removed by renormalization. A convenient way to implement the 
renormalization is to add counterterms to the action (|T0|): 



6S 



dt / d^x \ 5pi)^ - -5g{il)^f + ■ ■ 



(52) 



In perturbative calculations, the counterterms dp and 5g should be treated as quantum 
corrections. They can be expanded according to the order in the quantum correction: 



5p = 5ifi + S2P + ■ ■ 
Sg = hg + 52g + -- 



(53) 
(54) 



To obtain the free energy after the renormalizations of p and g, we substitute /i — > /i + (5/i 
and g —* g + Sg into ( |T7| ) and expand in the order of the quantum correction. The complete 
expression to second order in the quantum corrections is 



^(/i) 



9 



+ 



1 

4m 



(0,-1 



(4m/x) 



+ 



rngp ^ 2 
8 32^ -^' 



2-di/i + -irdig 



2/-i,-i/i,i - 3A%) + -/i,i(5i/i 



-2^52P+^52g 



\5,pf + 2^J,p5,g - ^{5,gf 



(55) 



9 9 9 9 9 

By repeating each of the steps in Section |3]^ including the effects of the counterterms, we 
obtain an expression for the energy density that takes into account the renormalization of p 
and g: 



£{p) 



^0 + -^9P^ + 



-—Io,^i{2mgp) 
Am 



C P T 

pdip+ —Oig 



+ 



mg'^p 
16 



9 



^+|2<^-'^ 



2L 



1,-1 



/i, 



il: 



p 



+ '-^h,iSig 



PS2P+ —S2g 



.(56) 
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The ultraviolet divergences in (^6|) that are independent of p can be cancelled by £q. The 
counterterms ^i/i, ^2/^, ^iQi and 829 can be determined by demanding that the p-dependent 
power ultraviolet divergences cancel. For example, the integral appearing in the first order 
quantum correction can be written 



-^0,-1 = / TTTT^ PyP"^ + 'imgp -p - mgp + 




+ ^9P - -^^j ■ (57) 

The first integral converges and each term in the second integral gives a power ultraviolet 
divergence. The integral of the p^ term in the second integral of (|57|) is independent of p 
and can be cancelled by ^(0). The remaining divergences can be cancelled in ( ^6|) by taking 

g r dPp 

*'» = —Iwjs?- '^^) 

For later convenience, we have generalized the integrals to an arbitrary number of spatial 
dimensions D. Similarly we can determine the counterterms 62P and 62g by demanding 
the cancellation of p-dependent power ultraviolet divergences in the second-order quantum 
corrections in (^): 

^2/^ = ^ ( / 7^^ 1 ( / 7^^:;. I , (60) 




^^9 = ^^ T^^i, ■ (61) 



In D = 3 dimensions, all the power ultraviolet divergences can be removed by renormal- 
izations of p, g, and Sq. However there is still a logarithmic ultraviolet divergence coming 
from the J term in (^61). That divergence is the only obstacle to completing our calculation 
of the energy density to second order in the quantum corrections. The predictive power 
of quantum field theory lies in the fact that the same renormalizations must remove the 
ultraviolet divergences from all physical quantities. The renormalizations that remove the 
ultraviolet divergences from the ground state energy density must also remove ultraviolet 
divergences from the amplitudes for the low-energy scattering of atoms in the vacuum. In the 
next Section, we calculate quantum corrections to the amplitudes for scattering of atoms in 
the vacuum. After identifying the renormalization that removes the logarithmic divergence 
from the energy density, we will complete the calculation of S in Section ^. 

4 Scattering of Atoms in the Vacuum 

In this Section, we calculate quantum corrections to the T-matrix elements for 2 — > 2 scat- 
tering and for 3^3 scattering of atoms in the vacuum. We determine the renormalizations 
that are necessary to remove ultraviolet divergences from these T-matrix elements. 

11 



4.1 Perturbative framework 

Atoms in the vacuum can be described by the action ([lOD with the chemical potential /i set 
to 0. A perturbative framework for calculating their scattering amplitudes can be obtained 
by separating the action into a free part and an interaction part as follows: 



5'free[^] 
SintH'] = I dt 



Jdt Jd^xtfj^Udt + ^jtlj, (62) 

' d'x ^-^g (V^t^)2 _ JL^3 (^t^)3 _ }j _ (63) 

We can read off the Feynman propagator from S'frce- Its Fourier transform is 

D(uj, k) = — ^ (64) 

This propagator is represented by a solid line with an arrow as illustrated in Fig. 0(a). The 
4-particle and 6-particle interactions in ( |63|) are represented by vertices connecting 4 lines 
and 6 lines, respectively, as shown in Fig. |^(b) and 0(c). These vertices have equal number 
of arrows entering and exiting. This reflects the conservation of the number of atoms, which 



follows from the phase symmetry ■?/' — > e*"'?/' of the action consisting of (|62D and (p3|) 



4.2 2 — > 2 scattering 

Two atoms with momenta ki and k2 can scatter into states with momenta k'^ and kg that 
are allowed by conservation of energy and momentum. The probability amplitude for the 
scattering process is given by the T-matrix element T(ki, k2; k'^, kg). 

The only terms in the action that contribute to 2 — >■ 2 scattering are those that are fourth 
order in ip. If the only such term is (V'^V')^) the scattering is purely S-wave. In this case, the 
T-matrix element is a function of a single variable: 

r(ki,k2;k;,k'2)=r(gi2), (65) 

where gi2 = |ki — k2|. The center-of-mass energy E is related to gi2 hy E = g^g/^^^i. 

The T-matrix element for 2 — >^ 2 scattering can be expanded in the order of quantum 
corrections: 

^(gi2) = -9 + T,{qu) + T2{qi2) + . . . . (66) 

The first term comes from the tree diagram in Fig. ^(a). The first quantum correction comes 
from the one- loop diagram in Fig. §(b). Using contour integration to evaluate the energy 
integral, we get 

TM2) = ^j (2,r)3p2-p.(ki+k2)+ki.k2-^6- ^^^^ 

The integral has a linear ultraviolet divergence, but the divergence can be removed by renor- 
malization of the coupling constant g. Including the counterterms from the tree diagram in 
Fig. ^(a), we get 

mg"^ f d^p 1 

Ti(gi2) = ^^ / yr-T^^ 7- — , . ^ , , -. --dig. (68) 

2 7 (27r)'^ p^ - p ■ (ki -f- k2) + ki ■ k2 - «e 

12 



To maintain rotational symmetry, we must shift the integration variable by p ^ p + (ki + 
k2)/2 before imposing an ultraviolet cutoff |p| < Aj/y. The resulting integral can be evalu- 
ated analytically and we obtain 

47r^ ibvr 

The counterterm 6ig, whose value was determined in (|59|), precisely cancels the linear ultra- 
violet divergence in (^). The final result is 

mg'^ 
^1(^12) = i-^qi2- (70) 

A particularly convenient method for regularizing ultraviolet divergent integrals is dimen- 
sional regularization. The number of spatial dimensions D is taken to be a complex variable. 
The integral is evaluated as a function of D in a region of the complex D-plane where it 
converges. This defines an analytic function of D which can be analytically continued to 
D = 3. After shifting the integration variable in (|67D by p ^ p + (ki + k2)/2 and then 
integrating over the angles in D dimensions, we obtain 

''''*^' = (4.)-/T(f)/o *p^-,y4-..- '''' 

The integral converges for ReZ) < 2 and is given by 

Analytically continuing to D = 3, we recover the result (|70D . 

One of the great advantages of dimensional regularization is that integrals that contain 
no scale are set identically equal to 0: 

'^''^ p" = 0. (73) 



(2vr) 



D 



This formula can be derived by first integrating over angles in D dimensions, and then 
separating the integral over p into two pieces corresponding to p < p* and p > p*. The 
integral over p < p* can be evaluated for D large enough that it is convergent in the in- 
frared. The integral over p > p* can be evaluated for D small enough that it is convergent 
in the ultraviolet. Upon analytically continuing the two integrals to D = 3, we find that 
they cancel exactly. Because of the identity ([73|), dimensional regularization sets pure power 
ultraviolet divergences to 0. Thus the counterterms 6fi and 6g, which are given by the inte- 
grals in (|55|)-(pTD, vanish. With dimensional regularization, the only ultraviolet divergences 
that require explicit renormalization are logarithmic divergences, which appear as poles in 
D — 3. Nontrivial counterterms are therefore needed only to cancel logarithmic ultraviolet 
divergences. 
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The quantum corrections to the scattering amphtude from higher order diagrams, hke 
the two-loop diagram in Fig. §(c), form a geometric series and can be summed up exactly. 
The nth term in the series is 

T„(gi2) = -g (-^^^12)" • (74) 

Summing up the geometric series, the complete 2 — > 2 scattering amplitude is 

^(^12) = , ^. ~\,,^ . . (75) 

1 + i'mgqi2/[l'o'K) 

The imaginary part of this T-matrix element is precisely that required by the optical theorem. 

4.3 3^3 Scattering 

Three atoms with momenta ki, k2 and ka can scatter into states with momenta k'^^, kg and 
kg that are allowed by conservation of energy and momentum. The probability amplitude 
for 3^3 scattering processes in which all 3 atoms participate is given by the connected 
T-matrix element, which we denote by T(ki, k2, ka; k'^, kg, kg). For simplicity, we consider 
only the center of momentum frame, where ki + k2 + ks = and we use the shorthand 

r(123 ^ 1'2'3') = T(ki, k2, ka; k^, k'^, k^) . (76) 

The connected T-matrix element for 3^3 scattering can be separated into the terms 
that involve a single virtual particle in the intermediate state and the remainder, which is 
called the one-particle-irreducible (IPI) part of T: 

T(123 -^ 1'2'3') = T^^\l23 -^ 1'2'3') 

777 

(123) (1^') ki ■ k2 - (ki + k2) • kg + fcg^ - te 

where gi2 = |ki — kg] and qi'2> = m — kg]. The sums are over cyclic permutations of ki, 
k2, and ka and of k'^, kg, and kg. The term in the sum that is given explicitly corresponds 
to the 2^2 scattering of particles 1 and 2 to produce particle 3' and a virtual particle. A 
subsequent 2^2 scattering of the virtual particle and particle 3 produces particles 1' and 
2'. Examples of diagrams that contribute to the sum are the tree diagram in Fig. ^(a) and 
the one-loop diagram in Fig. |^(b). 

The leading contributions to the IPI T-matrix element for 3^3 scattering come from 
one-loop diagrams like the one in Fig. |^(a). After using contour integration to integrate over 
the loop energy, we obtain 

r/P^(123 ^ 1'2'3') = -m'/ Y. E ^(123 ^ 1'2'3') , (78) 

(123) (1'2'3') 

where 

J(123 ^ 1'2'3') = J (2vr)3 (p2 + p . kg + k^ . k^ _ ,,)(p2 + p . ^ + k/ . k/ _ ,e) ' ^^^^ 
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This integral is ultraviolet convergent. 

The next most important quantum corrections to the IPI T-matrix element come from 
two- loop diagrams, such as those in Fig. ^b), ||(c), and ^(d), and from the insertion of a 
counterterm 6ig into the one- loop diagram in Fig. |](a). After using contour integration to 
integrate over the loop energies, we obtain 

r2iP^(123 ^ 1'2'3') = — r^'''\l23 -^ 1'2'3') 

+^VE E [Ti{qu) + T,{qr2')]I{l23^l'2'3') 

(123) {1'2'3') 

d^p f d?q 1 



+ -^VE E /■ 



(123) (1'2'3') 
X 



(27r)3 J (27r)3 p2 ^ g2 ^ r^ - 2mE - ie 



(p2 + p . ka + ki ■ k2 - ze)(g2 + q . k^ + k'l ■ k^ - ze) 

+ (p2 + p . ka + ki ■ k2 - It) (p2 + p . k^ + k'l ■ k^ - 2e) J ' ^^°^ 

where r = |p + q| and i? = {k\ + /c^ + k\)/2m is the total energy. The integral over q in 
the last term of ( [80| ) has a linear ultraviolet divergence. Using the expression (^) for the 
counterterm 5ig, we can see that the first term on the right side of ( ^0]) cancels the linear 
divergence from the integral over q in the last term. 

After cancellation of the linear divergence, there remains an overall logarithmic divergence 
in the integral over p and q in (|80D . This is evident from scaling p -^ tp and q -^ 
tq. As t ^ cxD, the integrand scales like 1/t^ while the integration volume scales like 
t^. The divergence is independent of the external momenta and therefore corresponds to 
a point interaction between the three particles. The divergence can be cancelled by the 
counterterm 5g-j, associated with the [ip'^ipY term in the action ([T0|) . However, if we include 
that counterterm, we must also for consistency include the contribution to T(123 -^ 1'2'3') 
from the coupling constant g^. Thus we must add to (|80|) the contribution from the tree 
diagram in Fig. ^(c): 

AT = -{g^ + 5g^). (81) 

We choose to use dimensional regularization to regularize the integral in (|80|). The 
logarithmic ultraviolet divergence then appears as a pole in Z) — 3. The pole in the integral 
over p and q in (^) is identical to that of the following integral, which is evaluated in the 
limit D — i> 3 in Appendix A. 4: 



d^p f d^q 



(277)^ J (27r)^ [ (p2 + g2 + ^2 + 2^2) (p2 + /,2)(g2 + ^2) 

1 
+ 



47r - 3^3 / 



(p2 + g2 + ^2 + 2fi;2)(p2 + ^2)2 2(g2 + n'^){p^ + ^2)2 j 

1.13459) fi:'(^-=^) . (82) 



1927r3 \D-3 



The last term in the integrand of ( P^ ) cancels the linear divergence in the integral over q of 
the previous term, but does not change the logarithmic ultraviolet divergence, which gives 
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the pole in D — 3. The pole term in (|80D is therefore 



T2^^\123 -^ 1'2'3') 



3(47r-3V3)^3 4 



pole 6An^{D - 3) 



m-'g'' . (83) 



The pole must be cancelled by the counterterm Sg^ in (|81|). One of the simplest renormaliza- 
tion schemes is minimal subtraction |jlO|. This scheme defines a running coupling constant 



g^{ti) that depends on an arbitrary renormalization scale k. The minimal subtraction pre- 
scription is to choose the counterterm to be a pure pole in D — 3 multiplied by a power of 
k: 

Sgsi^) = _ 3(4vr-3v/3) 4^,(^_3) 
^^^ ' 647r3(D-3) ^ ^ ^ 

The exponent of k is chosen so that both sides of (|4D have the same engineering dimensions 
even when D ^ 3. Without such a factor, renormalized quantities would involve logarithms 
of dimensionful quantities. The power of k in (^) is determined by dimensional analysis. 
With h set equal to 1, the terms in the action must be dimensionless. Time has dimensions 
mL'^ when ^ = 1, where L refers to length. The integration measure Jdtjd^x therefore 
has dimensions mL^^'^. Since ifj^^ip is a number density, if) has dimensions [^/^J = L^^/^. The 
dimensions of the coupling constants are then [g] = L^~'^/m and [g^] = L'^^^'^/m. The 
power of K in (34) provides the extra factor of L^"^^^^ required for dimensional consistency. 



Physical quantities cannot depend on the arbitrary parameter k introduced through the 
counterterm (|8^) . The coupling constant gsin) must therefore depend on k in such a way that 
the combination gs+Sgs is independent of k. This statement can be conveniently expressed in 
the form of a renormalization group equation. Using (|8^ , the condition {d/dti){g2, + 5g^) = 
reduces in the limit D ^ 3 to 

d 3(47r - 3^3) 3 4 

""T.'' = 32vr3 "^ ' ■ ^^^^ 

Since m and g are independent of k, the solution to the equation (^3]) is 



3(47r-3y3) 3 4, K' 



93{f^') = ^3(«:) + ^ ^_ 3^ ^ mV In - . (86) 

This equation tells us that the parameter g^ is a running coupling constant that varies 
logarithmically with the renormalization scale k. The renormalization scale k, can be inter- 
preted as the inverse of the spatial resolution. As k, increases, the spatial resolution becomes 
finer and part of the "pointlike" 3—^3 scattering amplitude represented by the diagram 
in Fig. ^(c) is resolved into the successive 2^2 scatterings represented by the diagrams 
in Figs. |^(c) and ^(d). The contributions from the two individual two-loop diagrams have 
opposite signs and the net effect is that the coupling constant g^ increases as k increases. 

The renormalized expression for 7^^^^(123 — > 1'2'3') is given by the limit as D — > 3 of 
the sum of (|80D and (p]). To take the limit, we must extract the pole in D — 3 from the 
integral over p and q in (pOD, so that it can be cancelled by the counterterm Sg^. This can 
be accomplished by subtracting the integrand on the left side of (|8^) from the last integrand 
in (pUD, and then adding to ( PPP the right side of ( P^ ) multiplied by 9m^g'^. The integral 
over p and q is then convergent in D = 3, but depends logarithmically on the scale k. 
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We have eliminated the logarithmic ultraviolet divergence from T^^^ by expressing it 
in terms of the renormalized coupling constant (73 (k) defined by dimensional regularization 
and minimal subtraction. The resulting expression for T^^^ can serve as a definition of 
gsin) that makes no reference to the regularization scheme. Someone who prefers a more 
physical definition of the coupling constant can define fi'3,phys to be equal to the value of 
— T^^^(123 —>■ 1'2'3') at their favorite configuration of the initial and final momenta. Any 
such coupling constant can be expressed in the form 

^3,phys = gsi^) + C{K)m^g\ (87) 

Since its definition makes no reference to the renormalization scale k, the physical coupling 
constant satisfies (d/dK.) g^^pi^ys = 0. The renormalization group equation (p5| ) then implies 
that the coefficient C(k) in ( ^7|) is a linear function of In^. It therefore vanishes for some 
value Kphys and we have 

5'3,phys = 5'3('^phys)- (88) 

Thus any physical definition of the coupling constant is equivalent to the running coupling 
constant gsin) evaluated at a particular value of the renormalization scale k. Thus there is 
little to be gained by using a more physical definition of the coupling constant. 

The integrals in 7^^^^ + AT are functions of the renormalization scale k, and the initial and 
final momenta. For momentum configurations in which the squares of the momenta and their 
inner products are all comparable in magnitude, the only scales in the integrands are k^ and 
mE, where E is the total center-of-mass energy. Since the integral varies logarithmically with 
K, it must also depend logarithmically on mE. If the ratio of these two scales is sufficiently 
large, the integral is dominated by the logarithm. However, since the dependence of the 
logarithm on k is cancelled by gsin), there must be a large cancelling contribution from 
the 5^3 term. The dominant terms in the second order quantum correction to the T-matrix 
element are therefore 

Tr + AT « -M.) - 5(11-|^™V In '^ . (89) 

The large logarithm can be avoided by choosing the renormalization scale of the running 
coupling constant (73 (k) to be k = ymE. Thus the most appropriate choice for the renor- 
malization scale K in the 3^3 scattering amplitude is the magnitude of the typical momenta 
of the scattering particles. This choice will avoid a large cancellation between the two terms 
infl). 



5 Renormalized Energy Density 

In this Section, we complete the calculation of the energy density in Section ^ by using the 
renormalization of gj, to remove the logarithmic ultraviolet divergences from the second-order 
quantum correction. Our final result is expressed in terms of parameters g and (73 (k) that 
can be defined purely in terms of the scattering of atoms. We then discuss the case of alkali 
atoms, and argue that in this case the dependence on (73 can be eliminated in favor of a 
logarithmic dependence on the length scale set by the van der Waals potential. 
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5.1 First-order quantum correction 

To simplify the calculation of the quantum corrections, we use dimensional regularization to 
regularize ultraviolet divergences and minimal subtraction to carry out the renormalization 
of the parameter g^. One of the great advantages of dimensional regularization is that it 
sets power ultraviolet divergence to 0. The only ultraviolet divergences that require explicit 
renormalization are logarithmic divergences, which appear as poles in D — 3. Since the 
ultraviolet divergences removed by the renormalization of /i and g are power divergences, the 
explicit renormalization of these parameters is unnecessary with dimensional regularization. 
This can be seen from the expressions (^8D-(0) for the counterterms, which all vanish by 
the identity (|73|). Thus the expression ( |^ ) for the energy density after renormalization of 
/x and g collapses to (|5T|) . The integration constant Sq, which is used to set £^(0) = 0, is 
also zero in dimensional regularization. This follows from the fact that the only scale in 
the dimensionally regularized integrals /„ „ and J is 2mgp. These integrals therefore vanish 
when p = 0, since there is no momentum scale. 

From (pTD, the first order quantum correction to the energy density is 

Slip) = ^10,-ii'^mgp) . (90) 

Using the expression for the dimensionally regularized integral given in ( |A12| ), this becomes 

SM^^I^M^. (91) 

Using (^ to express g in terms of the scattering length a, we recover the classic result for 
the first order quantum correction given in (p. 

5.2 Second-order quantum correction 

Dimensional regularization eliminates the power ultraviolet divergences from the second- 



order quantum correction in (|5lD , but the expression still contains a logarithmic ultraviolet 
divergence. The renormalization of this divergence requires including the effects of the 
coupling constant g^ and its counterterm at tree level. The term that must be added to the 



free energy density can be read off from the integrand in the action (|T0[): 

A^(/i) = hg, + 5gsW. (92) 

6b 

Following the effect of the correction through to the energy density, we find 

A^(p) = h93 + Sg3)p'. (93) 

6b 

With dimensional regularization, the complete second-order quantum correction is the 



sum of (p3) and the correction in (pl|): 



1 r / A , r / M 3 , ^9^P T , 9 



£2{p) = T^ [9.{^) + 5g,{n)] p' + -^J + ^ {l\_, - 2/_i,_,/i,i - II,) . (94) 



The values of the integrals J, /_i._i, and Ji i in dimensional regularization are given in ( |A2^ 
( |A13|) , and ( [A14| ) with A^ = 2mgp. The pole in D — 3 in J is cancelled by the pole in the 



counterterm Sg^ln) given in (Q), but the cancellation leaves a logarithm of 2mgp/ k^ in the 
limit D — > 3. Combining all the terms, our expression for £2 is 

SM = y.(^V + %^ (m ^ + O.80) ,n^oV . (95) 

where (73 (k) is the running coupling constant defined by the minimal subtraction renormaliza- 
tion prescription. The expression (|95D is independent of the arbitrary renormalization scale 
K. The renormalization group equation (jS^) implies that the explicit logarithmic dependence 
of (|95[) is cancelled by the ^-dependence of (73 (k). If h? differs by orders of magnitude from 
2mgp, there is a large cancellation between the logarithm in ( ^31) and the term containing 
(73 (k). Such a large cancellation can be avoided by choosing the renormalization scale to be 
K = y/2mgp. This is the momentum scale at which the dispersion relation for the Bogoliubov 
mode changes from linear to quadratic. Our final result for the energy density to second 
order in the quantum corrections is obtained by adding the corrections (|9ll) and (|95|) to the 
mean-field contribution: 

1_2 , 1 i2rngY^\,/2 



4 oOtt^ m 



1 
^ 36 



fe W + '*':: '^' f'..^ + 0.80 ),nV 



647r^ V K 



p\ (96) 



It would be very difficult to measure the coupling constant (73 (k) experimentally by 
studying the 3-body scattering of atoms. It would also be difficult to calculate this parameter 
theoretically from a microscopic description of the interaction between atoms. Thus (73 (k) 
must be treated as a phenomenological parameter. The predictive power of the result ( p6D 
for the energy density resides in the fact that the same coupling constant (73 (k) appears in 
the second order quantum corrections to other physical quantities, such as the dispersion 
relation for the Bogoliubov modes. All of the low-energy observables of the Rose gas can 
be calculated to second order in the quantum corrections in terms of two phenomenological 
parameters: the S-wave scattering length a and the coupling constant (73 (k) associated with 
3^3 scattering. 

5.3 Alkali atoms 

Our final result (^) for the energy density is of limited utility because it depends on the 
phenomenological parameter (73 (k). It may be possible to neglect the dependence on (73 if 
the second quantum correction is dominated by the logarithm. We argue that this is indeed 
the case for typical alkali atoms, and that the dependence on g^ can be eliminated in favor 
of a logarithmic dependence on the length scale set by the van der Waals interaction. 

The interaction between two atoms at low energy can be described by a potential V{r) 
that has a repulsive core at very short distances, an attractive region at short distances 
comparable to the size of the atom, and a long-range behavior given by the van der Waals 
potential: 

V{R) - -^ . (97) 

19 



The scattering length a is extremely sensitive to the short-distance behavior of the potential. 
Small variations in the depth or range of the potential can easily cause a to vary from — oo 
to +00. However, given a random distribution in one of the short-distance parameters of 
the potential, the distribution of a is concentrated in the region where \a\ is less than or 
comparable to the van der Waals length defined by |]11 



ma \^/^ , , 

^) ■ (^«' 

For most values of the short- distance parameters, a is comparable in magnitude to £y. 
However, a varies dramatically with the short- distance parameters near the critical values 
at which a new 2-body bound state appears. As a parameter passes through its critical 
value, a approaches ±cxd, changes discontinuously to =Foo, and then decreases in magnitude. 
The magnitude of a will be orders of magnitude larger than ly only if the short-distance 
parameter is in a narrow range around its critical value. We summarize this situation by 
saying that the natural magnitude of a is £y, and that much larger values require fine-tuning 
of the potential. 

One can use the concept of natural magnitudes to estimate the magnitude of the coupling 
constants in the lagrangian for an effective field theory. A coupling constant with dimensions 
L"- /m, where L refers to length, will by dimensional analysis have the form fiy/m, where 
/ is dimensionless. The assumption of naturalness is that the coefficient / is of order 1 
except when a short-distance parameter of the potential is tuned to within a narrow range 
of a critical value. Since the coupling constants g and gs in the action (p!OD have dimensions 
L/ni and L'^/m, respectively, the natural estimates for their magnitudes are \g\ ^ iy/^^ and 
I (73 1 ~ iy/rn, respectively. 

One can improve on these naive estimates by taking into account geometrical factors of 
4:71. For example, using the relation (||) between the coupling constant and the scattering 
length a, we obtain the estimate 

1 5' I natural ~ • 199; 

m 

We argue that the magnitude of g should be comparable to this natural estimate unless the 
2-body potential is tuned so that there is a bound state near threshold. 

We next consider the coupling constant gai^n) associated with 3 — > 3 scattering. This 
is a running coupling constant that depends on an arbitrary renormalization scale n. An 
estimate of the magnitude of this coupling constant must include a specification of the scale 
K, at which the estimate applies. Our estimates of natural values, which involve dimensional 
analysis, are based on the assumption that iy is the only important length scale. Thus the 
estimate must apply for momentum scales k that correspond to the length scale iy. We will 
therefore assume that the estimate of the natural value applies to 5^3 (k) for k comparable to 
271 /iy. Making a guess for the appropriate values of 47r, our estimate for the natural value 
is 

I93(.)l.,„, ~ ^ hr.^f. (100) 

m ly 

We will verify that this guess passes a simple consistency check. If both g and g^ have natural 
values, then we would expect their values not to change dramatically under changes of the 
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renorinalization scale by a factor of 2 or 3. Using the solution (|86[) to the renormalization 
group equation for (73 (k), we see that the change in g^ from the evolution of k by a factor of 

6 is 

3(47r-3V3) 04 ,^^^, 

A^3 ~ 327r3 ^^ ■ ^ °^ 

If (yf has the natural magnitude given in (^), the estimate (|100|) is approximately equal to 
( |101|) , which indicates that our guess of the factors of in in (|100|) is at least reasonable. 

We now consider an atom for which g^ has the natural magnitude given by ( P^OOj) . It 
could be much larger if the potentials describing 2-body and 3-body interactions are tuned 
so that there is a 3-body bound state near threshold. We cannot exclude such a possibility, 
but if we pick an alkali atom at random it is unlikely. We proceed to consider the three 
cases where g is much smaller than, comparable to, and much greater than the natural 
estimate given in (|99D . If g is much smaller than |5'|naturah then renormalization has little 
effect on the value of g^. The m^g'^ term in the second order correction to the energy 
density in (|96D is negligible compared to the g^ term. In this case, we cannot calculate the 
second order correction without knowing the value of g^. Next we consider the case where 
g is comparable to l^fl natural- The g^ term in (^) is then comparable in magnitude to the 
constant term multiplying m^g^. These two terms can be neglected only if the coherence 
length {2mgp)~^^'^ is orders of magnitude larger than £y/27r. In this case, the logarithmic 
term dominates and we obtain an estimate of the second order quantum correction that is 
independent of the unknown constant g^: 



Finally, we consider the case where g is much larger than If?! natural- In this case, (73 will 
quickly evolve under renormalization to a value comparable in magnitude to ( [101|) . Again 
we find that the g3 term in ( p6D is comparable in importance to the constant under the 
logarithm. If the coherence length is orders of magnitude larger than the van der Waals 
length, these terms are negligible compared to the logarithmic term and the expression ( p6D 
for the free energy density reduces again to ( |102| ). In summary, the approximation in 



is valid provided that g3 is not unnaturally large compared to the estimate ( p.OCI| ), that a is 
not unnaturally small compared to Cv-, and also that 

ln(pa4)| > 1. (103) 

Under these conditions, we can eliminate the dependence on 5^3 in favor of a logarithmic 
dependence on ly Using (^ to express ( |102| ) in terms of the scattering length a, we obtain 



the expression for the energy density given in (0). 

6 Conclusion 

We have calculated the second order quantum correction to the energy density of a homo- 
geneous Bose gas. This is the first correction that depends on an atomic physics parameter 
other than the S-wave scattering length a. We identify that parameter as a coupling con- 
stant g3 that specifies the point-like contribution to the 3^3 scattering of atoms in the 
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vacuum. The result for the energy density in terms of g and the running couphng constant 
gz{K) is given in (|96|). In the case of alkah atoms, we argued that the dependence on g^ can 
be ehminated in terms of a logarithmic dependence on the length scale Cy set by the van der 
Waals interaction. The resulting expression for the energy density is given in (H). 

Thus far the only alkali atoms for which Bose-Einstein condensation has been successfully 
carried out are ^^Rb, ^^Na, and ^Li. Our calculations apply only to atoms with positive 
scattering length, such as ^'''Rb and ^^Na. The scattering lengths for these atoms are 60±15 A 
for '^^Rb and 29 ±3 A for ^SNa |l|]. The parameter a in the van der Waals potential (|^ 
is roughly 7 keVA^ for Rb and 1.8 keVA^ for Na. The van der Waals length iy defined 
in (p8| ) is therefore about 60 A and 30 A, respectively. In both cases, the scattering length 



is comparable in magnitude to (.y, so that one of the conditions for (|102|) is satisfied. The 



condition ( p.03| ) depends on the density p and will be satisfied if p is orders of magnitude 
larger than l/(a£y). 

The magnitude of the quantum corrections increases with the number density. To see 
roughly how important these corrections are in existing magnetic traps, we evaluate them 
for the typical densities at the centers of the traps in the earliest experiments []T], |^ . For '^^Rb 
atoms with number density p = 3x 10^^/cm^, the correction factor in (H) is 1 + 0.004 — 0.0002. 
For ^^Na with p = 3 x 10^^/cm'^, the quantum correction factor in the energy density is 
1 + 0.01 — 0.002. In both cases, the second order quantum correction is an order of magnitude 
larger than one would have guessed by squaring the first-order quantum correction. The 
second order correction is relatively large because the logarithm in (0) is large, having the 
value -14 for ^'^Rb and -12 for ^SNa. Thus the condition ([T03|) for the validity of (D is 
indeed satisfied. 

The estimates given above suggest that the number densities in the Bose-Einstein con- 
densates that have been produced thus far are not sufficiently high for the effects of quantum 
fluctuations on the energy density to be measurable. Since the first quantum correction scales 
like y/p, the quantum corrections can be made larger by increasing the number of atoms in 
the trap. Unfortunately, the peak density p scales like N'^/^ [Q, so iV must be increased by 
orders of magnitude before the effects of quantum fiuctuations on the energy density will be 
measurable. There are however other observables that may be more sensitive to the effects 
of quantum fluctuations. These effects may also be more important at temperatures near 
the phase transition for Bose-Einstein condensation. We hope that our explicit calculation 
of second-order quantum corrections for the energy density of a homogeneous Bose gas will 
stimulate further work on quantifying the effects of quantum fluctuations on experimentally 
measurable observables. 
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A Loop Integrals for the Energy Density 

In this appendix, we evaluate the integrals that are needed to calculate the second order 
quantum corrections to the energy density. We also list the expressions for each of the 
two-loop diagrams. 

A.l Energy integrals 

The energy integrals can be evaluated by contour integration. The energy integral for the 
one-loop subdiagrams in Figs. ||(a)-H(c) are 

^^ 1 = L_ (M) 

271 uo^ - e\p) + le 2e{p) ' ^ ' 

where e{p) is the Bogoliubov dispersion relation given in (0). The energy integrals for the 
two-loop diagrams in Figs. |Kd)-|(g) are 

duJi f duj2 1 



27r J 27r [ul - e^{p) + ie] [u^ - e^{q) + ie] [{u^ + cua)^ - e^{r) + ie] 

^ (A2) 



4e{p)e{q)e{r) [e{p) + e{q) + e{r)] 
duo I f duj2 UJ1UJ2 



27r J 27r [cul - s^p) + ie][ujl - e\q) + ie][{uj^ + uj^f - e^{r) + ze] 

= CA3^ 

Ae{r) [e{p) + e{q) + e{r)] ' 

A. 2 One-loop momentum integrals 

The one-loop and two-loop corrections to the ground state energy density involve momentum 
integrals of the form 

where we have introduced the following notation for the integration over a momentum in D 
spatial dimensions: 

'"" (A5) 



P . (21,)" 



The integrals ( [A^ ) satisfy the identities 



d n 

;J-m,n = ~77 -'m+l,n+2 ) (A-O) 



c^A2 '"'" 2 

A J-m,n = -'m-l,n-2 ~ -'m+l,n • (,A7j 



After integrating over angles, the integral is 



1 /-oo lp2\m+(D-n-2)/2 

^m,n = (4^)D/2r(|) i, "^^^ (p2 + A2)n/2 ' (^8) 
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If the integral is convergent, we can use integration by parts to derive the identity 

{D + 2m- n)Im,n = nIm+2,n+2 ■ (A9) 



The integral ( [A8| ) is ultraviolet divergent ii D > 2{n — m) and infrared divergent if 
D < n — 2m. If the integral converges, its value is 

T p / 2n-2m-D N -p / D+2m~n \ 

J - ^ ^ 2 >^y 2 > KD+2{m-n) ( ^^V^\ 

'm,n - ^^^^^/2 r(f)r(f ) • ^ ' 

If the integral ( |A^ ) is ultraviolet or infrared divergent, it can be regularized using dimensional 
regularization. The regularized integral is obtained by analytically continuing the expression 
( |M0| ) to D = 3. The result is 



T -p / 2n—2m-3 \ p / 3+2m-n N 
T _ ^ ^y 2 )^y 2 I A 3+2(m~n) /AllA 

The one-loop integrals that appear in the ground state energy density are Iq-i, -?^-i,-i, and 
/i,i. In 3 dimensions, these integrals have power ultraviolet divergences. With dimensional 
regularization, they are given by the finite expressions 

^0,-1 = T^A', (A12) 

1-1,-1 = -Aa3, (A13) 

/i,i = ^A^ (AM) 

A. 3 Two-loop momentum integrals 

The two-loop correction to the ground state energy density (^TJ) involves a linear combination 
of the integrals 

J ^ r r {p/V^T7^)' (g/y/^TA^)" (r/y/f^+A^)" 

where r = |p-|-q|. In D = 3, these integrals have quartic and quadratic ultraviolet divergences 
that cancel in the combination of integrals J given in (^51). The expression for J can be 
written 



J ■ ■ 



P -'q pv^p^"+A^ + gVg2 + A^ + rVr^ + A^ 



6p 2 Vp2 + A2 



-\/p2 _|_ ^2 p 

.(A16) 



3pqr P\/q'^ + A^v^r^ + A^ 



Vp" + AVg' + AVr^ + A2 qry/p^ + K' 

This integral still has linear and logarithmic ultraviolet divergences. By subtracting and 
adding appropriate terms in the integrand of J, we can isolate the linear and logarithmic 
divergences into separate terms: 

J = J\in + <^log + <^num • (A17) 
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The term containing the hnear ultraviolet divergence is 



J,i 



lin 



2- 



Vp' + A' 



p 



qr 



(A18) 



The term in ( [A17 ) containing the logarithmic ultraviolet divergence is 

2 



Ji 



log 



-A^ 



p Jq [ (p2 + A2)(g2 + A2)(p2 + ^2 + ^2 + 2A2) 



+ 



p2 + g2 + ^2 + 2A2 2(g2 + A2) 



(p2 + A2)2 



(A19) 



The integral Jnum obtained by subtracting (|A18|) and (|A19|) from (|A16|) is convergent in 
D = 3 dimensions and can be evaluated numerically. It is convenient to symmetrize the 
integrand over the six permutations of p, q, and r in order to avoid cancellations between 
different regions of momentum space. The resulting expression is 



Jr\ 



p.q6 ^ 



1 



6p 



2Vp^T7^ 
(pqr) kPVp'^ + A2 + gVg2 + A2 + rVr^ + A^ [Vp^TA^ P 

3pqr p\/g^~+A^\/r2~+~/V2 



+ ^ 

gZ 



+ 



^p2 _^ A2^g2 _,_ A2v/r2 + A2 
p \/p^ + A2 



qr^/p^^^A^ 



Vp' + A2 



P 



2A^ 



(p2 + g2 + ^2 + 2A2)(p2 + A2) (g2 + A2) 



+ 



p2 + ^2 + ^2 + 2A2 2(g2 + A2) 



A^ 



(p2 + A2)2 



(A20) 



Since A is the only scale in the integrand, dimensional analysis implies that the integral is 
proportional to A^. Evaluating the coefficient of A^ numerically, we obtain 



Jn 



2.10 X 10"^ A^ 



(A21) 



Because of the severe cancellations between the various terms in the integral, we were only 
able to calculate it to 3 significant figures. 

We evaluate the ultraviolet divergent integrals Jun in ( |A18|) and Jiog in ( |A19|) using 
dimensional regularization. The integral over q in 
the integrand, and therefore Jiin 
the limit D 



T8|) vanishes since there is no scale in 
0. The integral ( |A19|) is evaluated in Appendix [A. 4 



m 



3, and the result is 



(A22) 



J,o, = ^^^^^ f^- 1.13459 V^+^^^-^^- 
^ 1927r3 \D-3 J 

Adding ( |A21| ) and (|A22D , we obtain the complete result for J using dimensional regulariza- 
tion: 



J = ^""-^f (^— + 0.57) A^+^(^-^) . 
1927r3 \D-3 J 



(A23) 
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A. 4 Evaluation of Jiog 

The term ( [A19| ), which contains the logarithmic ultraviolet divergence in the integral J, can 
be written 



'log 

where Ki and K2 are the following integrals: 
K2 



1 



P K (p2 + g2 + ^2 + 2A2) (p2 + A2) (^2 + ^2) 
1 1 



p-yq 



p2 + ^2 + ^2 + 2A2 2(g2 + A2 



(p2 + A2)2 



(A24) 

(A25) 
(A26) 



where r = |p + q|. 

We first consider the integral Ki. Setting r"^ = p^ + q^ + 2y> ■ f\ and then introducing 
Feynman parameters, the integral ( [A25| ) becomes 



Ki 



"^^io '^^ JpJc,[{l-y)p^ + {l-x)q^ + ZYi-ci + K'f 



(A27) 



where z = 1 — x — y. Dimensional regularization allows us to shift and rescale the momentum 
variables. We can eliminate the dot product from the denominator by making the shift 

p ^ p — 27]TrTCt- After rescaling p by (1 — y)^^^"^ and q by f *■ ~^-'^^~^^~^ / j ^ the integral 
factors into a Feynman parameter integral and an integral over the momenta: 



Ki 



dx 
Jo 



l~x 



dy[{l-x){l-y)-zy4] 



-D/2 



1 



/p7q(p2 + g2 + A2)3- 

The integral over the momenta can be evaluated analytically: 



(A28) 



p "'q 



(J92 + g2 + A2)3 



r{3-D] 



D 



^2{D-3) _ 



The gamma function has a pole at D = 3. To obtain Ki in the limit D 
expand the Feynman parameter integral in powers oi D — 3: 



(A29) 
3, we need to 



dx I dy 
Jo 



:i-x)il-y)-^- 



-3/2 



[l-x){l-y)~- 



-3/2 



In 



:i-x){l-y)-- 



(A30) 



where z = 1 — x — y. The first integral in (|A30|) can be computed analytically and has the 
value 47r/3. The second integral has to be computed numerically and has the value —9.43698. 
Extracting the pole in D — 3 from the gamma function in ( [A30| ) and keeping all terms that 
survive in the limit D -^ 3, we obtain 



Ki 



967r2 



^2(D-3) 



D-3 



+ 1.12646 + 7 -ln(47r) 



(A31) 
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where 7 is Euler's constant. 

We next consider the integral K2 in (|A26|) . By introducing a Feynman parameter, it can 
be written 



Ko 



dX (1 Xj / / I ~ — ~ ~ a~oTq 77 \ — o o a~oTq 1 ' y^oZi] 

Jp Jq \ [p^ + xq^ + a;p ■ q + A^J'^ [(1 — xjp'^ + xq'^ + A^\-^ j 



By shifting and rescahng the momentum variables, we can reduce the integral over the 
momenta to ([A29D . In the first term of ( |A32| ) , we shift p ^ p — f q and then rescale q 

by \ ^ ) • In the second term, we rescale p by (1 — x)^^/^ and q by x^^^"^. After 
integrating over p and q, we obtain 



K, 






2(47r 



XjX 



-D/2 



(1 - x/4)-^/2 - (1 - x)-^/^ 



(A33) 



To obtain K2 in the limit D —^ 3, we need to expand the integrand in ( |A33|) in powers of 
D-3: 



K, 



n^-^l^^io-^ir,. 



2(47r; 



D 



:i-x) 



1 



D-3 



dx 



)(l-x) x(4-x) 



^x^{A-xf y/x3(l -x) 
1 



a;^(4 — xY 



4 



x^fl — x) 



ln[x(l — x)] 



(A34) 



The integrals can be evaluated analytically. Extracting the pole in D — 3 from the gamma 
function and keeping all terms that survive in the limit D — >^ 3, we obtain 



K. 



647r3 



^2(D-3) 



1 



47r 



1 3 

_ -In 1 + 7- ln(47r) 

D-3 3^3 2 4 ^ ^ ' 



(A35) 



A. 5 Two-loop diagrams for the Free Energy Density 

The two-loop vacuum diagrams that contribute to — 'ji72(/i, vq) are shown in Fig. | 
diagrams can be reduced to momentum integrals by integrating over the loop ener^ 
the identities in Appendix A.l . Expressed in terms of the integrals Im,n and Ji^m, 
in (|A4j) and ( |A15|) , the contributions of the individual diagrams to r22(/x, wo) are 



fi(") = 


3 r2 


Qf = 


= ^^9l-i,-ih,i 


n? = 


3 ,2 

= 64^'-'- 


nf - 


J 2—2 T 


a-) - 


3 2—2 T 



|. These 
ies using 
, defined 



(A36) 
(A37) 
(A38) 
(A39) 
(A40) 
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^F^ = -^rng'vlJ^^,_,^^ (A41) 

Qi'^ = -lmg'vlJ.,.0,0. (A42) 

o 

Adding up these diagrams, we get (pl). 
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Figure Captions 

1. Propagators and interaction vertices for the real- valued fields ^ and rj. 

2. Two-loop vacuum diagrams that contribute to the thermodynamical potential. 

3. Propagator and interaction vertices for the complex-valued field ip. 

4. Diagrams for 2—^2 scattering: the tree- level diagram (a), a one-loop diagram (b), and 
a two-loop diagram (c). 

5. One-particle- reducible diagrams for 3^3 scattering: a tree- level diagram (a) and a 
one-loop diagram (b). 

6. One-particle- irreducible (IPI) diagrams for 3^3 scattering: a one- loop diagram (a) 
and two- loop diagrams (b), (c), and (d). 
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